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^^^ ' The dynamics of the 1 + 1 abehan Higgs model with fcrmions is studied in the large 

Cf^ , Nf approximation, on a real-time lattice. The Bose fields obey effective classical 

equations of motion which include the fermion back reaction. The dynamics of the 
quantized fermion field is treated with a mode function expansion. Numerical re- 
sults are shown for renormalizability, nonequilibrium dynamics and the anomalous 
charge, and Pauli blocking. 
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PsJ ' 1 Introduction 

o' 

0^ ' Non-perturbative approximations to real-time dynamics in Quantum Field 

^^ I Theory include large TV, Hartree or semiclassical approximationsm Fields are 

written as a sum of a mean or classical field and a quantum part (treated in 

OhI practice with a mode function expansion), and effective equations of motion 

f-^ • can be derived that couple them. Such an approximation includes e.g. the 

g_) \ Hard Thermal Loops in 3 -I- 1 dimensional hot gauge theories. Furthermore, 

T^ • the effective equations contain the divergences of the quantum theory, and 

these can be renormalized in the usual way. 

However, in actual numerical implementations, the emphasis has been on 
, ; , homogeneous mean ficldsH This excludes dynamics related to inhomogeneous 

C^ ' classical configurations such as kinks, sphalerons, and bubbles. Furthermore, 

when the quantum mode functions interact via a homogeneous mean field, 
there is no 'mode mixing'. This makes the long-time behaviour of the system 
peculiar. Allowing for inhomogeneous mean fields immediately removes the last 
two objections. The drawback is that the numerical problem becomes much 
more demanding. Space dependence of the mode functions is no longer given 
by plane waves, but is determined by the full partial differential equations. 
Here we'll discuss results for a i -f 1 dimensional model, where the numerics 
can be handled with brute force. A moie detailed account, and much more 
technical details, can be found elsewhereH 
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2 Abelian Higgs model with fermions and effective equations of 
motion 

A toy model for electroweak baryogenesis is the 1 + 1 dimensional abelian Higgs 
model with axially coupled fermions. The action is given by 



S=- d'x 



As the electroweak theory, this model has a non-trivial bosonic vacuum, spha- 
lerons, and an anomalous global symmetry: fermion number violation. 

We want to study the dynamics numerically in real-time. We find it then 
convenient to use a real-time lattice formulation, for several reasons: the lattice 
acts as a gauge invariant ultraviolet cutoff, there are exact symmetries for finite 
lattice spacing, and the resulting integration algorithm is simple and stable. 
Fermions on a lattice give of course rise to the fermion doubler problem. To 
deal with this, we use Wilson's fermion method in space, and interpret the 
doubler in time as a second flavour. For technical reasons, it is convenient 
to transform the model to one with a vector gauge symmetry, by performing 
charge conjugation on the right-handed fermions only. The anomalous current 
is then axial, with the anomaly equation Qz{t) — Q^iQ) = C{t) — C(0), where 
C = — X^x^i/^"" ^^ ^^^ Chern-Simons number. Finally, it is also convenient 
to use a real 4-component Majorana field ^. 

Effective equations of motion can be derived from firsi- principles, by du- 
plicating the fermion field Nf times, and taking Nf -^ ooBfl In the language 
of the Introduction, the Bose fields represent the mean fields, and the fermions 
the fluctuations. The resulting bosonic equations are (in the temporal gauge 
Ao=0) 
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These are similar to the classical equations, but with the fermion back reac- 
tion present, i.e. the fermion current {j'i) ~ j{'^'^(3j'^P2'^) and force (F) = 
j{^(3{pi + ips)^) (in continuum notation), pi are matrices appearing in the 
Majorana description. The fermions are treated with a mode function expan- 
sion, which in the continuum would read 





Figure 1: Chern-Simons number C for three values of the cutoff, A'^ = 16, 32, 48 and G/e = 
0.5. Left: the bare v'^ = 10 is kept fixed. Right: the renormalized v^ = 11.25 is kept fixed. 



a labels the complete set of mode functions (it typically contains momentum), 
the creation/annihilation operators are time independent and determine the 
initial quantum state (we use vacuum, {h]j)a) = 0), and the mode functions 
themselves are solutions of the Dirac equation 



'^i{do + d'^)U^{x,i) =nD{A,<i>)U^{x,t), 



(3) 



with Tio{A, (j>) the Dirac hamiltonian in presence of the Bose fields. Note that 
we now have a closed set of equations. For a precise treatment of the lattice 
fermion doublers m time, and the initial conditions for the mode functions, we 
refer to our paperH The initial conditions for the Bose fields are such that only 
long wave lengths are excited, and that Gauss' law (lb) is satisfied. 



3 Renormalizability, nonequilibrium dynamics, and Pauli blocking 

The equations contain the ultraviolet divergences of the quantum theory, and 
need to be renormalized. In the full quantum theory only the scalar self energy 
is divergent, and because we only integrated out the fermions, we expect to 
find only the fermion loop contribution. Let's take a closer look to the scalar 
field equation (||). If we take for simplicity (j) — vr/^/?, (i.e. equal to the 
renormalized v.e.v., vb is the bare parameter), and Ai = 0, and also evaluate 
the fermions in this background, (0) reduces to the gap equation 



X(vl-vl)^^G{F)ivj,X {F){vn) = ^ (^^ logiV + O(iVo) 



(4) 
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Figure 2: Left: Energy of the Bose fields, fermion field, and the conserved sum. Right: 
Chern-Sinions number (dashed), and the anomalous charge (solid), for G/e = 0. 



An explicit calculation of {F) shows that it is indeed a logarithmically divergent 
sum, indicated with log A^ in <^p) {N is the number of spatial lattice points). 
This divergence is canceled by the appropriate vg. In practice, we fix Vfj, and 
then find (for certain N, i.e. lattice spacing) the corresponding bare parameter 
vb, from (^). The result of this procedure is demonstrated in fig. |l|. As 
expected, a proper renormalization gives converging physical results. 

In fig. |l| the Chern-Simons number undergoes a damped oscillation, with a 
small amplitude. Remember that the sphaleron sits at half-integer C. To probe 
this non-perturbative physics, we put more energy in the initial Bose fields. 
The results are shown in fig. y. The initial state is clearly out-of-equilibrium 
and there is energy transfer from the Bose fields to the fermions. The Chern- 
Simons number is initially very wild, but as the bosonic energy decreases, the 
Bose fields start to feel the sphaleron barriers, and plateaus become visible. 
The anomalous charge follows the Chern-Simons number, in accordance with 
the anomaly equation: AQs = AC. 

Energy transfer from the Bose fields to the fermionic degrees of freedom 
cannot go on forever because of Pauli blocking. This becomes clear if we choose 
the initial energy in the Bose fields to be much larger than before (compare 
figs. H, H). Again there is rapid energy transfer (until approx. et = 1000), 
but then there is a rather sharp transition after which the energy in both 
subsystems remains approximately constant. A heuristic picture is as follows: 
the state with maximal fermion energy is given by a completely filled Dirac 
sea and sky. Note that this situation is not physical, it is a manifestation of 
the fact that for finite N we only have a finite number of states. Nevertheless, 
a calculation of the 'free' fermion energy (i.e. neglecting the fluctuating Bose 
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Figure 3: Left: as in fig. Q (left), with G/e = 0.1. Right: blowup of the fcrmion energy, the 
'free' maximal fermion energy is indicated with a dashed line. 



fields, keeping only — V]i/\/2) when all states in the spectrum are completely 
filled, gives for our choice of parameters and for A^ = 32, E^^^ /e = 815. 
Indeed, a blowup shows that the fermion energy fluctuates around this value. 



4 Summary 

As an approximation to non-perturbative dynamics in quantum field theory, we 
studied a coupled system of classical Bose fields and a quantized fermion field. 
Skipping many technical details, we presented three numerical results. Many 
more questions can be addressed in this framework (e.g. long time behaviour, 
inclusion of CP violation, sphaleron rate in the presence of fermions), and we 
hope to report on these issues in the future. 
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